Meta-nematic, smectic and crystalline phases of dipolar fermions in an optical lattice 
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It has been suggested that some strongly correlated matter might be understood qualitatively in 
terms of liquid crystalline phases intervening between the Fermi gas and the Wigner crystal or Mott 
insulator. We propose a tunable realisation of this soft quantum matter physics in an ultra-cold gas. 
It uses optical lattices and dipolar interactions to realise a particularly simple model. Our analysis 
reveals a rich phase diagram featuring a meta-nematic transition where the Fermi liquid changes 
dimensionality; a smectic phase (stripes); and a crystalline, 'checkerboard' phase. 



A picture of strong correlations [l[ unfolds as follows: 
as the strength of correlation increases, the Fermi gas 
condenses into a more correlated Fermi liquid. At this 
phase transition, the Fermi surface may change volume 
or even 0, [H alter its topology. Then, the first 'electronic 
liquid crystal' state forms: the nematic Fermi liquid 
accessed through a symmetry-breaking deformation of 
the Fermi surface (a Pomeranchuk instability (5(). As the 
strength of correlation increases further, a smectic phase 
develops. In this 'striped' phase the Fermi liquid state 
is lost as electrons localise, but only in one direction. In 
the other direction the stripes behave as Luttinger liq- 
uids. Thus such liquid crystalline phases are intimately 
related to dimensional crossover phenomena (a subject 
of intense current interest, both in a condensed matter 
context |6| and in cold atomic gases |7J). Eventually, in 
the limit of very strong interactions, the particles localise 
completely, forming a Wigner crystal or a Mott insulator. 

Experimental evidence for abrupt changes of Fermi 
surface volume or topology exists for heavy fermions Q. 
A nematic state is supported by transport measurements 
in YBa2Cu30g+y 0] (with the transition rounded by lat- 
tice anisotropy). There is evidence of nematic order in 
quantum Hall devices [lj| ■ A Pomeranchuk instability 
may explain 'hidden' order in the he avy fermion URu2Si2 
111 ] and the ruthenate Sr 3 Ru2 07 [12j. Smectic phases 
exist in manganites [T3| and cuprates [3]. In summary, 
there is evidence that elements of the scenario in Refs. m 
resemble the physics of strong correlations. Yet in order 
to establish its general usefulness, a system that can be 
tuned from the Fermi gas all the way to the localised state 
and is amenable to theoretical treatment is necessary. 

In recent years it has become possible to realise strong 
correlations in highly-tunable cold atom experiments 
[HI]. Simple models, like the Hubbard model, can be 
realised precisely. Unfortunately even the 2D Hubbard 
model is very difficult to solve, even approximately. Ex- 
periments of that type must therefore be regarded as 
'quantum analogue simulations' fl6j | . Here we propose an 
optical lattice set-up featuring dipolar fermions in an ex- 
ternal field. The system consists of a 2D stack of chains, 
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Figure I: (color online) Right: Proposed experimental set- 
up. Dipolar atoms or molecules (1-6) are loaded on a 2D, 
anisotropic optical lattice in a strong external magnetic or 
electric field (B). The field is oriented so that all interactions 
between lattice sites are repulsive and there are no intra-chain 
interactions (see text). Left: Interaction potential in recipro- 
cal space for three values of the anisotropy parameter a — 0.5 
(top), 1.0 and 2.0 (bottom). 



each of them containing free fermions. In the absence of 
interactions the ground state is a non-interacting Fermi 
gas with a nearly flat Fermi surface. Using an external 
field to produce a particular orientation of the dipoles 
relative to the lattice, we introduce a strictly inter-chain 
interaction as a perturbation, and address the stability 
of the ID Fermi surface with respect to it. We argue that 
the system will feature a meta-nematic transition where 
the quasi-lD Fermi surface becomes fully 2D, competing 
with phase transitions into smectic and crystalline order. 

The combination of optical lattices with dipolar inter- 
actions can be used to realise exotic Hamiltonians with 
novel states. This has been discussed extensively for 
bosons 17, 18, 2(| 21, [22|. Here we propose to use 
a fermionic isotope with a large magnetic dipole mo- 
ment such as 53 Cr [23| . Alternatively, dipolar molecules 
[2I . I22I or atoms cooled into a Rydberg state [25| may 
be used. Dipolar interactions between fermions in an ex- 
ternal field are expected to display a range of interesting 
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phenomena, including Fermi surface deformations (in 3D ^ 



traps) 26], exotic quantum Hall states (in rotating 2D 
traps) [27j and a 'super Tonks-Girardeu regime' (in ID 
traps, for either bosons or fermions) 28]. 

Proposed experimental setup. — We propose to com- 
bine a polarizing external field (electric or magnetic, 
depending on whether we are exploiting the magnetic 
dipole of atoms or the electric dipole of molecules, for 
example) with an anisotropic, 2D optical lattice, as il- 
lustrated in Fig. [TJ For sufficiently intense lasers, the 
lattice is in the tight-binding limit with one orbital per 
site [l5|. By allowing three different intensities for the 
three pairs of lasers we can make the system completely 
2D and create "chains" along which hopping can occur, 
while keeping hopping perpendicular to the chains much 
smaller. Bonds along these two directions are represented 
by the thick and thin lines in Fig. respectively. Tun- 
ing the effective wavelengths of the two in-plane lasers 
provides independent control of the lattice anisotropy 
a = a^/aj_. Isotropic scaling of the lattice controls the 
relative strength of interactions. 

We shall suppress the interaction between atoms or 
molecules that are on the same chain (for the sake of 
brevity, in what follows we assume the case of atoms 
in a magnetic field without loss of generality). To this 
end we exploit the dependence of the dipole-dipole in- 
teraction, V{R) = d 2 [l - 3cos 2 6»]/ |R| 3 , on the angle 
9 between the vector giving the relative positions of the 
two dipoles, R, and an external field, B strong enough 
to fully polarise all the atoms. We have represented this 
depedence schematically in Fig. [1] by the line L, disc D 
and cone C around atom number 1. In these directions 
the interaction with another atom is maximally attra- 
tive [9 = 0], maximally repulsive [9 = tt/2] and null 
[6 = arccos (l/\/3) ~ 54.736°], respectively. In the pro- 
posed arrangement (see figure) the applied field is at pre- 
cisely this latter 'magic angle' to the chains and perpen- 
dicular to the inter-chain bond direction. Thus atom 1 
does not interact with other atoms on the same chain, 
such as 2 and 4. Since on-site interactions are forbid- 
den by Pauli's exclusion principle, in this setup there 
are no intra-chain interactions [391 ] . Moreover inter-chain 
interactions are always repulsive and maximum in the 
direction of perpendicular hopping. The strongest repul- 
sion corresponds to the closest sites on the two adjacent 
chains: V (R) = d 2 /a\ (e.g. atom 3). Interactions with 
other sites on the two adjacent chains (e.g. 6) can also 
be made comparatively weak (see below). 

Model. — The single-particle Hamiltonian is H^ ^ = 
~Yh,i (*Fi,z^+i,( + t±4,A,l+i +H.c.) , where ct ; cre- 
ates a fermion on the i^ 1 site of the chain, t\i is 
the intra-chaing hopping amplitude and t± <C tu is 
the inter-chain hopping. Defining the operator creat- 
ing a fermion with wave vector k = (kn,k±) by — 
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we obtain H^ Qp = J2k e k £ l £ k, 
which has an almost flat Fermi surface given by the ze- 
roes of the 'bare' dispersion relation, = — 2tn cos (fcii) — 
2tj_ cos (kj_) — ji (/x is the chemical potential). 

The dipolar interaction has to be evaluated at the lat- 
tice sites. In the configuration discussed above, for rela- 
tive coordinates x = a\\i and y = a±l (in units of the lat- 
tice constans in the parallel and perpendicular directions) 

it gives Vij cc y 2 / [x 2 + y 2 ] . Fig. Q] shows the Fourier 
transform of this interaction, which depends strongly on 
the anisotropy ratio a. For a > 2 it is well approx- 
imated by V(kn,k±) sw 2V cos(fci) i.e. nearest neigh- 
bour only interaction. First we restrict ourselves to this 
limit. The interaction part of the Hamiltonian is thus 
Ant = ^i,M,i+i^i,l+i^i,h and the full form is 



^-^hop + ^int' 



(1) 



The three parameters controlling our model are ji/tn, 
t±/tn, and V/tn . Below we discuss possible ground states. 

Meta-nematic phase transition. — We start by evalu- 
ating the stability of the Fermi surface shape. We use as 
a trial ground state a Slater determinant of plane waves, 
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Hartree-Fock mean field theory is similar to those used 
to study Pomeranchuk 0, [3(| and topological 0] Fermi 
surface shape instabilities. 

The momentum distribution corresponds to a non- 
interacting Fermi gas with a renormalizcd dispersion re- 
lation e£ = —2t|i cos (fey) — 2t* ± cos (k±) — fi*. The struc- 
ture of -ffj n t is such that only the perpendicular hopping 
changes. It is given by the self-consistency equation 
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Numerical solutions are shown in Fig. [5] As the bare 
inter-chain hopping t± is increased, its renormalised 
value, t* ± , initially increases linearly but then has two 
bifurcation points, between which lies a first-order jump. 
The corresponding phase diagram is shown in Fig. [3] (note 
that for Fig. [5] we chose a very large value of V/tn, for 
clarity; for smaller values, the results are qualitatively 
the same, but the jump of t\_ is much smaller). The 
order parameter of this phase transition is the amount 
of delocalisation in the perpendicular direction, tp = 
(cJjC^z+i + H.c). We refer to the jump of tp as we 
vary t± as a meta-nematic transition in analogy with 
meta-magnetism (where the magnetisation jumps under 
an applied magnetic field). As V — > 0, the meta-nematic 
transition becomes more and more weakly first order and 
requires a larger value of t ± . At V = there is no longer 
a first-order transition, but the phenomenon survives at 
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Figure 2: (color online) Dependence of the renormalised 
transverse hopping t± on its bare value t± for V = 3tii and 
/i = — 1.5iii. The solid (dashed) lines correspond to solutions 
to the self-consistency equation @ that minimize (maximize) 
the energy. Insets: Fermi surface (a) just to the left of the 
bifurcation region and (b) just to the right, as indicated; and 
(c) dependence of the solutions on V for t± = (rightmost 
curve), 0.05t|| and 0.1t\\ (leftmost curve). 
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Figure 3: (color online) (a) Phase diagram of the Hamil- 
tonian of Eq. ([l]) for /i = ±1.9tii. The circles track the two 
bifurcation points ('a' and 'b' in Fig. [2} of the first-order meta- 
nematic quantum phase transition between the quasi- ID and 
2D phases. The solid line marks the line of quantum critical 
points separating these Fermi liquid states from the crystalline 
state (C). (b) Critical value of the lattice anisotropy a\\/a± 
for the dominant instability to be towards the crystalline (C) 
or smectic phases (S) as a function of the relative strength of 
the renormalised perpendicular hopping, t* ± , for /i* = ±1.9ti| . 



t± = t\\ +fi/2 as a 'two-and-half order Lifshitz transition 
[3l| (while remaining frist order for any V > 0). 

This quasi-lD to 2D transition induced by inter-chain 
interactions is in some sense the opposite of confinement 
[32! (a quasi-lD to ID transition induced by mtra-chain 
interactions [13]). A similar phenomenon is believed to 
occur in stacks of integer quantum Hall systems [34j |. 
where the chiral Luttinger liquids on the edges couple 
together, creating a 2D Fermi surface (the chiral Fermi 
liquid). In our cold-atom set-up, the meta-nematic tran- 
sition results from enhanced scattering when the poten- 
tial reaches the singularities at the edges of the ID bands. 
This is a density of states effect and hence we expect it 
to be robust to quantum fluctuations present for large 
values of V/tii and not taken into account by our mean 
field theory. It could be induced by changing the inten- 
sity of one of the lasers to tune t±, and detected by direct 
imaging of the Fermi surface [35| . 

Crystallisation. — The meta-nematic transition is not 
the only one possible in the sysem described by Eq. |T]). 
Interchain backscattering can lead to 'charge de nsity 
wave' (CDW) instabilities at low temperatures [361 ]. 
We probe the potential CDW instability by examin- 
ing the Fourier transform of the dynamic susceptibil- 
ity, X(k,t) = i($ |l>(k,t)pt(k,0)| Here p(k) = 
E q c q c q-k is the Fourier transform of the local occupa- 
tion number, and p(k, t) is its Heisenberg representation. 
The 'noninteracting' susceptibility (Lindhard function) 

is given by X (k,w) = / J!, q ~+«. +k in terms of the 

v ' q q+k 

renormalised dispersion relation e£. Treating the inter- 
action within the random phase approximation (RPA) 
gives X(k,u) = X (k, w)/[l + 2V cos(ki_)X {k, J)]. An 
instability at wave vector k occurs if the static compo- 



nent diverges, X(ki\,k±,u) = 0) — * 00. 

For t* ± -c tn, Xo(k, uj) is strongly peaked at (2kF,Tr) 
due to the strong nesting of the quasi-lD Fermi surface. 
Thus the system is unstable to a CDW of that periodicity 
at a critical coupling V given the Stoner criterion 1 = 
2VX (2k F ,n,Lu = 0). In the limit t* ± 0, the Fermi 
surface is perfectly nested and the peak in X becomes a 
logarithmic divergence, implying V — > 0. More generally 
one has to evaluate A"o(2fc^,7r,w = 0) to obtain V via 
the Stoner criterion. The results are plotted in Fig[3](a). 

Smectic phase. — Within the nearest-neighbor approx- 
imation for the interaction we have been using so far, 
the Fourier transform of the interaction potential V(k) 
is independent of kii , hence the dominant CDW instabil- 
ity is always at the peak of the Lindhard function, i.e. 
fen — 2kp. However, we now consider the full structure 
of the dipolc interaction. In particular, as seen in Fig. [T] 
when the ratio of lattice spacings a is not large, V(k) 
acquires a large dependence on kn . Hence so long as 
t* ± 7^ (i.e. -Xo(fc|| = 2k f) is finite), then there is a level 
of anisotropy of the lattice where the leading instability 
is at (0,7r) and not (2/cf,7t). 

The (0, 7r) instability is still a form of CDW - however 
as it breaks lattice symmetry in one direction only, it has 
smectic order. Figj3] (b) shows which instability takes 
place first as the overall strength of the interaction in- 
creases. Smectic order is favored when the fermions can 
lower their energy by crowding every other chain, paying 
a penalty in kinetic energy but taking advantage of the 
absence of intra-chain interactions. Note that the strong- 
coupling limit ground state is always a density wave. 

In the limit t± — > 0, the system becomes ID and one 
can employ bosonization. The bosonized Hamiltonian 
features backscattering terms responsible for the crys- 
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tallisation at arbitrarily small coupling [36(. If we artif- 
ically turn off these backscattering terms we find a sec- 
ond instability characterised by a softening of the holon 
dispersion relation. The critical coupling for this sec- 
ond instability coincides with the divergence of the RPA 
susceptibility at (0, n). Such accuracy of RPA is a re- 
sult of a special feature of Hamiltonian |T]), namely that 
interactions and single-particle dispersion take place in 
perpendicular directions in the limit t± — > 0. Indeed as 
a result the lowest-order vertex corrections vanish iden- 
tically in that limit. As tj_ grows such corrections will 
become increasingly important. For tj_ ~ fn strong cor- 
relation effects will modify the phase diagram, at least 
quantititatively. Describing these effects is beyond the 
scope of the present calculation. 

In summary, we have described a way to combine dipo- 
lar fermions in an external field with an optical lattice to 
realise a model featuring an array of chains with strictly 
inter-chain interactions. We have shown that the model 
has very rich physics, featuring competition between itin- 
eracy and localisation and between quasi one- and 2D be- 
haviour. The possibility to realise meta-nematic, smec- 
tic, and crystallisation transitions in a regime where the 
system can be described using Hartree-Fock and RPA, to- 
gether with the ability to introduce stronger correlations 
continuosuly into the system (for example, by tuning tj_ 
towards larger values) makes the proposed experiment a 
testbed for scenarios of correlated behaviour. 
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